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The nature of mobility of ions and water molecules in dilute aque- 
ous solutions of electrolytes (at most fifteen water molecules per 
ion) is investigated. It is shown that the behavior of the mobility 
coefficients of water molecules and ions, as well as the self-diffusion 
coefficients of water molecules, are determined by the radii of their 
hard shells rather than by the effect of the hydrogen bond net- 
work. It is established that the influence of hydration effects on 
the density of the system and the self-diffusion coefficients of water 
molecules does not exceed several per cent. Based on microscopic 
concepts, it is shown that the different behaviors of a K+ cation 
and an F~ anion with equal rigid radii are in good agreement with 
specific features of the intermolecular interaction described by the 
generalized Stillinger-David potential 1,2. 



1. Introduction 

For today, a substantial body of information is accumu- 
lated on the thermal motion of water molecules and ions 
in water electrolyte solutions (see [51-fTB]). Now, of inter- 
est is the experimental study of the kinetic and electro- 
physical properties of electrolytes. The specific features 
of motion of water molecules and ions have a certain 
effect on the thermodynamic properties of liquids and 
liquid systems [2] . The study of a structure and a char- 
acter of thermal molecular motion in electrolytes with 
the help of the quasielastic scattering of slow neutrons 

Table 1. Self-diffusion coefficients of water molecules 
in water solutions of singly charged electrolytes [5J 









LiBr(30) 


Lil(24.8) 


Di™' X 10^, cm^/s 2.3 2.18 






NaCl(15.9) 


NaBr(16.5) 


Nal(16.1) 


d)^^ X 10^, cm2/s 2.14 2.26 2.38 




KF(15) 


KC1(16.1) 


KBr(16.1) 


KI(16.4) 


Di™' X 10^, cm^/s 1.99 2.44 2.68 2.8 



[5l410| and molecular dynamics methods |11H13] acquires 
still larger importance. 

There also exists a large number of substantial theo- 
retical works [15,16, devoted to the description of prop- 
erties of electrolytes. The main efforts in them are con- 
centrated on the adequate allowance for hydration ef- 
fects and the correct reflection of the role of hydrogen 
bonds. However, a lot of questions (very simple, at first 
sight) remain practically unanswered. A typical example 
is peculiarities of the manifestation of hydration effects 
depending on the radii of ions and the signs of their 
charges. 

Let us illustrate this fact by the example concern- 
ing the behavior of the self-diffusion coefficients of wa- 
ter molecules as functions of the cation and anion sizes. 
The corresponding values of the self-diffusion coefficients 
Ds of water molecules in several dilute solutions of 
electrolytes at the temperature T ~ 296 K are given in 
Table 1 presenting the dependence of Dg on the anion 
size. The solution concentration in the table is speci- 
fied by the number Zw of water molecules per one ion 
(the value of z^ is given in parentheses near the chem- 
ical formula of an electrolyte, for example NaCl(15.9)). 
The minimal value z^ = 15 thus corresponds to the elec- 
trolyte concentration equal to 3.3 mole mutual effect of 
cations and anions can be neglected to a satisfactory 
accuracy. It is also worth adding to the table the self- 
diffusion coefficient of water molecules in the Csl(17.4) 

solution: i:)i'"^=2.71xl05 cm^/s. 

We note that, according to |17|[TS] . the self-diffusion 
coefficient of molecules in water at the same temperature 
T = 296 K is equal to 

D^'"'> = 2.35 X 10"^-- 



^cnr/s. 



(1) 



The number of hydrogen bonds formed by one 
molecule in water at the same temperature is approx- 
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imately equal to nn ~ 2.5 [T^H^ . In this case, the 
average distance between water molecules is close to the 
hydrogen bond length In « 2.8 A. Comparing M with 
the values of the self-diffusion coefficients given in Table 

1, one can conclude that the self-diffusion coefficients of 
water molecules considerably depend on the combina- 
tions of cations and anions even in dilute solutions of 
electrolytes. 

In this work, we analyze the mobility coefficients of 
water molecules as functions of the cation and anion sizes 
in the close connection with the behavior of the solution 
density. In this way, we try to make certain conclusions 
concerning the role of hydration effects in solutions of 
electrolytes. In addition, to discover the nature of asym- 
metry effects in the interaction of anions and cations 
with water molecules, the microscopic ideas about the 
character of their interaction are involved. This aspect 
of the problem is solved on the basis of the generalized 
potential of interparticle interaction proposed in ^In2p23| . 

2. Dependence of the Self-Diffusion and 
Mobility Coefficients of Water Molecules and 
the Mobility Coefficients of Ions on Their 

Size 

To analyze the results presented in Table 1, let us con- 
sider how the self-diffusion coefficients of water molecules 
in solutions of electrolytes correlate with the radii of dis- 
solved ions. The latter are not determined unambigu- 
ously, but considerably depend on the method of their 
determination. In this connection, we consider the fol- 
lowing ionic radii: 1) those found from crystallographic 
conceptions; 2) those chosen to favor the correct repro- 
duction of the molecular dynamics in solutions of elec- 
trolytes by means of computer simulation, and 3) those 
estimated from the ion polarization. To make the picture 
complete, we consider also the ionic radii determined 
from hydrodynamic considerations. The corresponding 
results are gathered in Table 2. 

In the first row of Table 2, one can see the crystallo- 
graphic radii Tc [24 . The second row contains the ionic 
radii r^ found from computer experiments |25| for the 

Table 2. Rigid (three upper rows) and Stokes radii 
of cations and anions 





Li+ 


Na+ 


K+ 


Cs+ 


F- 


c\- 


Br- 


J- 


re, A 


0.6 


0.95 


1.33 


1.69 


1.36 


1.81 


1.95 


2.16 


r,T, A 


0.76 


1.3 


1.67 


1.94 


1.56 


2.2 


2.27 


2.59 


Ta, A 


0.45 


1.12 


1.41 


2.02 


1.51 


2.33 


2.55 


2.93 


ri''\ A 


2.38 


1.84 


1.25 


1.19 


1.66 


1.21 


1.18 


1.19 


ri''\ A 


1.91 


1.88 


1.14 


1.15 


1.77 


1.33 


1.22 


1.35 



description of the dispersion (van der Waals) interaction 
between ions and water molecules. The ionic radii Tq, 
determined from the ion polarizations a (see Table 3) 
according to the formula 



= 1.5a 



1/3 



(2) 



are presented in the third row. It is worth noting that, in 
the simplest model of an ion presenting it as a conduct- 
ing sphere, the polarization is equal to a = ''sphj where 
Tsph is the sphere radius. With regard for the fact that 
a conducting sphere reflects ion properties only approxi- 
mately, the coefficient 1.5 in Q is chosen to provide the 
minimum deviation of the ionic radii r^ from r^- The 
radii re, r^, and r^ will be called rigid ionic radii. 

For the sake of comparison, the fourth and fifth rows 
of Table 2 present the Stokes ionic radii determined from 
the formulas 



rM 



and 



67r77/j,(-f) 



(3) 



riD) _ 



fceT 



Q'K'qD 



ii) 



,W 



where fi^^> and Dg^' stand for the mobility and the 
self-diffusion coefficient of ions determined by computer 
modeling in |25fl27| , whereas rj is the shear dynamic vis- 
cosity of electrolyte solutions ^28J . 

An important supplement to these data can be the 
Stokes radii of water molecules in electrolyte solutions 
gathered in Table 1 . They are determined by the formula 



rM 



knT 



GnriD, 



{w) 



(4) 



In contrast to r, and ri . 



and form Table 4 as well, 
the values of Ts practically do not depend on the type 
of a dissolved electrolyte, as is worth expecting. At the 
same time, the obtained radii rg are noticeably smaller 
than the hard-sphere radius of water molecules, if identi- 
fying the latter with a half of the hydrogen bond length 
[29. ,30] . Here, it is worth noting that hydrodynamic 
concepts can be used only on the scales significantly ex- 
ceeding the corresponding molecular ones. That is why 

Table 3. Ion polarizations, A^ |26| 



Li+ 


Na+ 


K+ 


Cs+ 


F- 


cr 


Br- 


I- 


0.029 


0.179 


0.83 


2.42 


1.04 


3.66 


4.77 


7.1 
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we make conclusion that the appHcation of hydrody- 
namic conceptions to the self-difFusion process of water 
molecules is inconsistent. 

One can see that the ionic radii used in the computer 
experiment somewhat exceed their crystallographic val- 
ues. In the both cases, one observes, however, a 
monotonous growth of the cation and anion radii with 
increase in masses. The same dependence is registered 
for the ionic radii obtained from the ion polarizations. 

On the contrary, the Stokes radii of water ions and 
molecules demonstrate the opposite behavior. The 
largest Stokes radius corresponds to a small Li+ cation, 
whereas the least one is observed for an I~ anion. There 
also exists a certain asymmetry in the behavior of cations 
and anions of the same radius. Indeed, the differences in 
the Stokes radii of K+ and F~ ions, as well as in their 
mobilities (see Table 5), significantly exceed the differ- 
ences in the values of their rigid radii r^, Va, and Tq,. 
This fact is to a certain extent unexpected. It should be 
related to the noticeably different character of the inter- 
action of water molecules with K"*" and F~ ions, as well 
as with all other cations and anions. This important 
fact is thoroughly analyzed in Section 4 by the example 
of K+ andF-. 

As follows from Tables 1 and 2, the self-diffusion co- 
efficients of water molecules as functions of their rigid 
radii satisfy the following regularities: 

1) for dilute lithium and sodium solutions of elec- 
trolytes, in which Tc < 2^' ^^^ following inequality 

is met: Z?s (el) < Ds ■ The condition Tc < ^Ir is 
evidently violated only for Cs"*". In solutions of potas- 
sium electrolytes, where rc(K"'") ^ ^/h, one observes a 
transition from the previous relation between the self- 
diffusion coefficients of water molecules to the inequality 
Di^\el)>Di^^; 

2) in solutions of electrolytes with a fixed cation, ex- 
cept for lithium ones, the self-diffusion coefficients of wa- 
ter molecules grow with the anion radius; 

3) in lithium electrolytes, the character of the depen- 
dence of Ds' (el) on the anion radius is opposite to the 
second conclusion. 

With regard for the fact that the concentrations of 
different singly charged cations and anions are close to 
each other, the differences in the behavior of Dg '(el) 

Table 4. Stokes radii of water molecules in solutions 
of electrolytes gathered in Table 1 

LiBr I Lil I NaCl I Nal I KF I KCl I KBr I KI 



can be caused by the following reasons: 1) geometri- 
cal factors and 2) their different effect on the structure 
of the local environment (hydration effects). The first 
possibility should be discarded, as the geometrical ob- 
stacles must decrease with the cation radius, which con- 
flicts with experimental data. But a reconfiguration of 
the local structure of water in the direct environment 
of cations and anions appears more noticeable, as their 
rigid radii increase. 



Indeed, comparing the values of r, 



M JD) 



and r. 



(w) 



one can see that, as the rigid radii of cations and anions 
increase, their Stokes radii approach the Stokes radius 
of water molecules. It is natural to interpret this fact as 
the joint drift of large ions together with water molecules 
that enter hydrated shells formed around them. It is 
worth noting that the effect of a local reconfiguration of 
the solution structure is insignificant, because the rise 
or fall of the self-diffusion coefficient of water molecules 
usually does not exceed ten per cent and is proportional 
to the mole concentration of electrolyte impurities. 

3. Analysis of the Hydration Energy of Ions 
and the Density of Electrolyte Solutions 

The dependence of the Stokes ionic radii on the rigid 
sizes of cations and anions is closely related to the phe- 
nomenon of the so-called positive or negative ion hydra- 
tion |14 | [3T | l32|. Indeed, according to |33j, the sign of 
the ion hydration energy is determined by the quantity 
AW — W — W\, where W is the average binding energy 
of two neighboring water molecules, and Wi is the aver- 
age value of their interaction energy in the presence of 
ions. An ion close to the neighboring water molecules 
results in their additional polarization and an increase 
of the attractive forces acting between them. In the ab- 
sence of hydrogen bonds, the difference W — Wi must 
be positive, as the both terms W and Wi are nega- 
tive. If a cation or an anion partially destroys hydro- 
gen bonds, then the absolute value of Wi can become 
lower than that of W, i.e. the negative sign of AW 
will testify to the local breakage of the hydrogen bond 

Table 5. Mobility coefficients of ions determined 
from the conductivity of electrolytes (/^c = /ic X 10^) | 27| 
and by computer simulation of the mean-square shift of 
an ion (/xr ^ Ar X 10^) and its autocorrelation velocity 
function (/x„ = /x„ X 10^) 125| 



J™) 



1.009 0.97 0.997 0.957 1.053 0.936 0.917 0.912 
0.934 1.025 1.015 0.951 1.034 0.949 0.881 0.85 





Li+ 


Na+ 


K+ 


Cs+ 


F- 


cr 


Br- 


I- 


Ac 


4.01 


5.19 


7.62 




5.7 


7.91 


8.13 


7.96 


Ar 


4.75 


4.98 


7.12 


7.32 


4.04 


6.88 


7.2 


6.23 


fl-v 


4.59 


5.02 


7.20 


7.36 


3.85 


6.42 


6.85 


6.27 
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network. Table 6 presents the dimensionless combina- 
tion AW = AW/3kBT,n (where 3A:bT„i is the energy 
of translational thermal motion of two water molecules 
at the crystallization temperature Tm — 273 K) as a 
function of the type of ions. Analyzing Table 6, one 
can see that Li+ and Na+ cations, whose rigid radii are 
less than a half of the hydrogen bond length (i.e. a 
half of the average intermolecular distance between wa- 
ter molecules) , conserve a rather intact structure of the 
local hydrogen bond network and intensify the interac- 
tion between neighboring molecules. A relatively small 
value of AW(Na+) can be interpreted as a consequence 
of the more significant bending of hydrogen bonds as 
compared to that as a reaction to the introduction of a 
small-sized Li+ cation. A K+ cation, whose rigid radius 
is close to a half of the hydrogen bond length, causes a 
damage of the local hydrogen bond network, which re- 
sults in the negative sign of AW. In the same way. one 
can explain the negative sign of AW for a Cl^ anion. 

At the same time, it is worth noting a significant asym- 
metry in the behavior of AW for K+ and F^ ions having 
almost identical rigid radii (see Table 2). The positive 
sign of AVl^(F~) testifies to the fact that an F^ an- 
ion does not completely break the local configuration 
of hydrogen bonds, though significantly deforms it. It is 
clearly proved by a relatively small value of AW^. To our 
opinion, the appearance of this asymmetry is caused by 
the polarization component of the interaction between 
ions and water molecules. The existence of the asymme- 
try in the behavior of cations and anions was first noted 
in [31IS5]. 

On the other hand, the role of hydration effects can 
be estimated, by analyzing the density of weak solutions 
of electrolytes (see Table 7) . We must actually estimate 
a deviation of the density of the real solution from that 
of the ideal one. By definition, an electrolyte solution 
can be considered as ideal if the volumes occupied in it 
by water molecules and dissolved ions acquire the same 
values as in water and electrolyte melts. For such an 
electrolyte solution, the following relation must be met: 



(0) 



vl">n^ + iviy'+vnnei 



,(0)] 



1. 



(5) 



Here, n^ and Uei denote the number densities of wa- 
ter and electrolyte molecules in the real solution, while 



,(0) ,,(0) 



,(0) 



Vw' , Vc ' , and vi are the volumes occupied by water 



Table 


6. Dependence of A.W on the type of ions [33j 




Li+ 


Na+ 


K+ 


F- 


cr 


AW 


3.4 0.06 -4.6 0.52 -2.0 



molecules and ions in water and in melted or solid elec- 
trolytes. The values of n^ and rid can be obtained from 
the definition of the mass density of an electrolyte solu- 
tion: 

p = mw^w + (mc -I- ma)nei 
and its weight concentration 



(TOc + ma)Wcl 

m^n^ + (toc + ma)nci 



Hence, 



(1-x)- 



rici = X- 



771c + "TlsL 



(6) 



A real electrolyte solution is described by the relation 
similar to (Isl: 



Ww»^w + (Wc + Va)nci = 1, 



(7) 



though it contains the real volumes occupied by a water 
molecule, cation, and anion. That is why a deviation of 
the dimensionless combination 



<5 = t;io)nw + (w(°)+t;i°))nei-l 



(8) 



from zero serves as the measure of adeviation of a water 
electrolyte solution from the ideal one. Moreover, pos- 
itive values of the non-ideality parameter (S > 0) will 
correspond to the formation of hydrated shells around 
ions, in which the water density will be larger as com- 
pared to the density of bulk water. Indeed, the formation 
of a hydrated shell is accompanied by a decrease of the 
volume occupied by a water molecule (u^ < Vw ). The 
variation of the volumes occupied by cations and anions 
in weak electrolyte solutions can be neglected. 

Let us substitute the following values of the parame- 



ters Uw •• , Vc , 



and Va in (18 1 



Pel 



where p- 



(0) 



,(0) 



rriw 

(0): 



,(") I „,(") 



and p'^^ 



p. 



are the water and elec- 



trolyte densities in the liquid and amorphous states. 
With regard for the number densities of water and elec- 
trolyte impurities determined by relations m\, Eq. (pi) 

Table 7. Densities of water electrolyte solutions at 
the fixed concentration x^i ^ 4 wt. 



cr 



Br- 



Li+ 

Na+ 
K+ 
Cs+ 



-0.001 
0.003 

0.005 



0.003 



0.004 
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takes the form 



P 



P_ 

(0) 



1 






1 



(9) 



The non-ideahty parameter of a dilute solution can be 
interpreted in somewhat another way. Taking into ac- 
count that the dominant contribution to 5 is given by a 
variation of the volume occupied by one water molecule, 
Eqs. (|8| and (|7| yield: 6 ~ n^Sv^. So, to an acceptable 
accuracy, the non-ideality parameter for a dilute solution 
is 

The densities of some electrolyte solutions, as well as the 
parameters of their non-ideality 5, are given in Tables 7 
and 8. It is worth noting that the non-ideality parame- 
ters obtained according to formula (|9| are averaged over 
the number of water molecules surrounding a certain ion. 
The change of the relative volume of water molecules in 
the first hydrated layer of an ion is larger by a factor of 
-^w/^i, where zi is the number of water molecules getting 
to the first hydrated layer. 

In particular, the non- ideality parameter for a NaCl 
solution at Xd = 4 wt. reaches 5 ~ 0.003. One 
can see that the relatively small values of mobility of 
lithium cations and self-diffusion coefficients of water 
molecules in lithium electrolytes correspond to negative 
non-ideality parameters. This testifies to the fact that 
lithium cations do not favor the formation of hydrated 
shells around them with densities exceeding that of wa- 
ter. In other cases, such hydrated shells are formed. At 
the same time, the smallness of 5 is an evidence of the 
fact that hydration effects in dilute electrolyte solutions 
have a weak effect on their densities. 

At relatively small concentrations of electrolytes, the 
deviation of a solution from ideality is proportional to 
its mole concentration c = ^ ^'"^ x: 



fa ec. 

The quantity s also can be interpreted as a measure of 
non-ideality of a solution. For the above-considered so- 
lutions, it lies in the range 0.1 < e < 0.5. 

Table 8. Non-ideality parameters of the solutions at 
the fixed concentration x^i = 4 wt. 





F- 


cr 


Br- 


I- 


Li+ 




-0.001 






Na+ 




0.003 






K+ 






0.003 


0.004 


Cs+ 




0.005 







The result described in two previous sections allow 
one to conclude that lithium cations do not radically 
change the local structure of surrounding water formed 
with the help of hydrogen bonds. As a result, their 
thermal motion can be considered as a drift in tempo- 
rary "cavities", whose generation is favored by the lo- 
cal hydrogen bond network. As the characteristic time 
of such drift, one can consider the time of existence 
of hexagonal rings in water, which should be identi- 
fied, according to [30], with the duration of the set- 
tled life of water molecules tq. The characteristic shift 
of a lithium cation can be estimated by the formula 
|Ar(Li+)| ^ ^6D,{U+)to, or |Ar(Li+)| « VT(Li+)To. 
The value of tq at the temperature T = 295 K is approxi- 
mately equal to (0.8-^1) x IQ-^^ s [TF. The self-diffusion 
coefficient of a Li+ cation is estimated with the use of the 
formula Ds{Li'^) — k-oTjlc x 10^. As a result, the charac- 
teristic shift |Ar(Li+)| of a lithium cation during an ele- 
mentary diffusion act is equal to |Ar(Li+)| sa 2.8 x 10^^ 
cm. The same order of magnitude is characteristic of 
the combination 'L'7'(Li^)To. This allows one to conclude 
that |Ar(Li+)| practically coincides with the average dis- 
tance between neighboring water molecules or the size of 
a hexagonal ring. A similar conclusion can be also valid 
as regards the thermal motion of a Na+ cation. Indeed, 
|Ar(Na+)|/|Ar(Li+)| ^ VAc(Na+)//ic(Li+) ^ 1.1. 

On the other hand, cations and anions, whose rigid 
radii exceed a half of the hydrogen bond length (n > 
I Zh), destroy the hydrogen bond network around them- 
selves and thus move freer. The rigid radii of K"^ and 
F^ ions are approximately equal to a half of the hydro- 
gen bond length (n « ^ Zh). That is why the features of 
their thermal motion considerably depend on the pecu- 
liarities of ion-molecule interaction. 

4. Interpretation of the Obtained Results Based 
on Microscopic Concepts 

To provide a clearer interpretation of the obtained re- 
sults, let us apply the microscopic approach, in which 
the interaction of water molecules and ions is described 
on the basis of the generalized Stillinger-David polariza- 
tion potential (2l[23] having the structure 



"SwI = $1 + $11 + $ 



III 



$ 



IV, 



(11) 



where $i, i = I, II, III, IV are the components of the po- 
tential that describe: $1 - direct Coulomb interaction 
with oxygen and hydrogens of a water molecule; <I>ii - 
the charge-dipole interaction between 1) the ion and oxy- 
gen and 2) the ion and hydrogens of a water molecule 
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RKrf = 4,0 A 




Fig. 1. Molecular system of coordinates 

(it is assumed that oxygen and the ion acquire a dipole 
moment due to their polarization); $ni - interaction of 
the dipole moments of oxygen and the ion, and $iv - 
repulsion between electron shells of oxygen and the ion. 
The explicit form of the contributions $i - $iv is given 
in Appendix. 

The positions of ions on equidistant surfaces is deter- 
mined by the polar (9) and azimuth {(/)) angles specified 
in the molecular system of coordinates (MSC) (Fig. 1). 
It is assumed that oxygen and hydrogens are located in 
the plane {x,y) of the MSC so that the axis y coincides 
with the bisector of the angle formed by O^^ and H^, 
HJ" and is directed opposite to hydrogens. 

The angular dependence of the energy $^1 ('i'wi = 
'i'wi/ k^Tm, where T^ — 273 K is the melting tempera- 
ture) of interaction between a water molecule and K+, 
F~ ions with almost coinciding rigid radii (rc(K+) = 
1.33 A, ra(F~) — 1.36 A), is calculated with the help 
of formula (fTTl and depicted in Fig. 2. Considerably dif- 
ferent anisotropy characters of the equidistant surfaces of 
the interaction of K+ and F^ ions with a water molecule 
are determined by the different structures of their elec- 
tron shells: for K+. the shell is argon-like, whereas it is 
neon-like for F^. This is the reason for the noticeably 
different polarizations of K+ and F~ ions (see Table 3) . 

At the same time, the situation somewhat simplifies 
if ions get to water. This is related to the fact that the 
interaction of an ion with surrounding water molecules 
is not accompanied by significant changes of the local 
arrangement of water molecules. Indeed, the average 
distance between molecules of water in a vicinity of its 
crystallization point is close to 3 A, whereas the sum 
of the rigid radii of water molecules is approximately 
equal to the hydrogen bond length, Zh ~ 2.8 A. So, the 
volume of "empty space" in water is insignificant. The 
introduction of K+ and F^ ions is actually possible if 
they substitute one of the water molecules in local con- 



figurations formed, first of all, by hard cores of the latter. 
It is worth noting that such cores are close to spheres to 
an acceptable accuracy [291150]. Hydrogen bonds give 
rise to considerable dipole correlations, as well as mul- 
tipole correlations of higher order. Moreover, structural 
violations in water surroundings will grow with increas- 
ing deviations of the cation and anion radii from those 
of K+ and F~ . Such structural perturbations of water 
cannot be interpreted, if one is based only on concepts of 
the dominant role of the hydrogen bond network. How- 
ever, the analysis based on the averaged intermolecular 
interaction between an ion and water molecules is quite 
sound. The interaction is naturally averaged due to the 
thermal rotation of water molecules. By definition, the 
averaged interaction potential Gwi {r) is described by the 
formula 



Gwi(0 



/ / $wi (r, 6*, ,?!))e-''*-i sin 9 d9 d<p 
//e-'3*wisin6'd6l# ' 



(12) 



where /3 — l/k-gT. The averaged interaction energy as 
a function of the distance between K+ and F~ ions and 
water oxygen is presented in Fig. 3, where we also show 
the asymptotic behavior of the averaged potential of in- 
teraction between an ion and a water molecule (dotted 
line) that can be calculated in the explicit form. 

Indeed, at large distances between an ion and oxygen, 
the energy of their interaction is determined by the main 
terms of the multipole expansion 



$wi(r) 



gildw -r) 



, (dw •r)(di -r) 



(13) 



which includes the charge-dipole and dipole-dipole in- 
teractions. Let the dipole moment of a water molecule 
be fixed and change only in direction, while let the ion 
dipole moment arise due to its polarization by the elec- 
tric field of the water dipole: 



oi 



, (dw ■r)r 



(14) 



Thus, the main contributions to the energy of dipole-ion 
interaction take the form: 



$wi(r) 



91 (dv 






r)2 



(15) 



At sufficiently large distances between a water molecule 
and an ion, 

exp(-^$wi) = l-/3$wi + --- ■ 
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K^ 



H 



i© 





a b 

Fig. 2. Equidistant surfaces (roi = 4.0 A) of the dimensionless energy of interaction #„! between a water molecule and K+ (a), F~ (b) 



That is why the asymptotic behavior of Gwi('") will be 
determined by the contributions 



Gwi(r) = ($wi)o-^($wi)o + 



(16) 



where the index "0" denotes the averaging over the 
isotropic orientation distribution of water molecules. 



Combining ( 14 1 . ( 15 1 , and ( 16 ) , it is easy to make sure 



that the main contributions to the averaged potential are 
determined by the following terms: 



Gwi(r) = -13 



2^2 



qfd. 



,aidv 



(17) 



One can see that, regardless of the sign of the ion 
charge, the averaged interaction between ions and a wa- 
ter molecule at large distances is of attractive character. 
For singly charged ions, the main asymptotic contribu- 
tion (from the charge-dipole interaction) does not de- 
pend on the type of ion. At the same time, the second 
contribution is not universal: it directly depends on the 
polarization of ions. 

Particularly, Eq. (17 1 together with the analysis of 



Table 3 make one expecting that, at rather large dis- 
tances (r > 10 A) between an ion and a water molecule, 
GwF-('') > GwK+('')- This conclusion agrees with the 
calculation results. At the same time, in the case of 
relatively small distances (r < 10 A), the averaged po- 
tential of interaction between an F^ anion and a water 
molecule appears deeper than Gwk+('') (see Fig. 3), ow- 
ing to which an F~ anion is less mobile than a K+ cation, 
which also completely agrees with the data in Table 5. 



20 


!3 -20 

O -40 

-60 



-100 




Fig. 3. Distance dependence of the averaged potential G„i{r) of 
interaction of a water molecule with K+ (a = 0.83 A'^) and F^ 
{a = 1.04 A3) 

5. Discussion of the Results 

The study of the transfer processes in dilute solutions 
of electrolytes performed in this work allows us to make 
the following conclusions: 1) the key role in the forma- 
tion of kinetic properties of water electrolyte solutions 
(first of all, in the behavior of the mobility coefficients 
of ions and water molecules) is played by their hard-core 
radii; 2) the standard understanding of the ion motion 
in "cavities" formed by the hydrogen bond network is, to 
our mind, inconsistent; 3) the hydrogen bond network 
in water is not determinative for the character of the ion 
drift motion, as the sizes of all cations and anions ex- 
ceed those of "cavities" of the local water structure. In 
addition, it should be taken into account that a certain 
local configuration of water molecules exists only during 
the time coinciding with the duration of the settled life 
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of molecules tq by order of magnitude. Speaking about 
the role of hydrogen bonds in the self-diffusion problem, 
we actually mean the manifestation of strong orienta- 
tion correlations between the dipole moments of water 
molecules and multipole moments of higher orders; 4) 
the behavior of the mobility coefficients of cations and 
anions in water agrees with that of the averaged poten- 
tials of interaction between ions and water molecules. 
Though this fact is established only for one cation and 
one anion, it still can be generalized to all other ions 
based on the similarity principle. 

Analyzing the density of dilute solutions of elec- 
trolytes, one can see that, if one ion is added per approx- 
imately 15 water molecules, the change of the relative 
volume occupied by one water molecule does not exceed 
a half of one per cent. This testifies to the fact that hy- 
dration effects are related, first of all, to the variation of 
the local orientation ordering of water molecules rather 
than to a local change of the density of the system. 

The established connection between the behavior of 
cations and anions and the behavior of the averaged po- 
tential is a result of the existence of thermal rotational 
motion of water molecules. This fact is most obviously 
manifested in the temperature dependences of the di- 
electric relaxation time and the shear viscosity of water. 



of the polarized ion I. The polarization of oxygen results in the 
appearance of the dipole moment do that characterizes the degree 
of deformation of the oxygen electron shells. The ion is polarized 
under the action of the hydrogen field of charges and the charge of 
oxygen in the water molecule. The polarization of the ion results 
in the appearance of the dipole moment dj characterizing the de- 
gree of deformation of the ion electron shells. With regard for the 
polarization contribution, one obtains 



-I-T 



(do ■ roi)gi 



''OI 



.7 = 1,2 



(di • ry )<jj 



[1 - L{roi)] + 



ll-K{nj)], 



(di ■ rio)<Jo 



^3 



[1 - L(no)]+ 



{D2) 



where do and di are the dipole moments of oxygen of the water 
molecule and the polarized ion I, respectively, whereas 1 — L{r) and 
1 — K{r) are screening functions (see [Tl[2ll23| V The dipole moment 
do is determined in the molecular system of coordinates, whose 
origin coincides with the center of mass of oxygen of the water 
molecule. The dipole moment dj is determined in the molecular 
system of coordinates, whose origin coincides with the center of 
mass of the ion. At large distances, the polarization contribution 
'I'ni is reduced to the interaction between the polarized oxygen 
of the water molecule and the ion charge, as well as between the 
polarized ion and the charges of the water molecule. 

The component ^m describing the dipole-dipole interaction of 
the ion and oxygen has the form 



#111 = <J>d(do, di)[l - K{roj/a)], 



(D3) 



APPENDIX 

Let us consider, in brief, the structure of each contribution to the 
energy of interaction between an ion and a water molecule (Eq. 
(n\). The first term $i P|2l[23] determines the direct Coulomb 
interaction between an ion (hereinafter denoted by character I) 
and oxygen and hydrogens of a water molecule: 



i=2 



*T 



E 



91 gj 

ri 



(DI) 



=0 'J 



where qi is the charge of ion I, j = 0,1,2, the index j numbers 
the charges of oxygen and hydrogen of the water molecule {j = 
corresponds to the charge of oxygen, j = 1,2 - to those of 
hydrogens). It is assumed that the charges are measured in units 
of the electron charge and gj = ±n {n is the ion charge), go = ~2, 
and qi = q2 = 1- 

The second contribution <l?n determines the potential of inter- 
action between the point ion charge and the polarized oxygen of 
the water molecule. Oxygen is polarized under the action of the 
field of charges of hydrogens in the molecule itself and the charge 

Table 9. Parameters of the potential [2l l32| of inter- 
action between a water molecule and an ion 





ki 


P3, A-i 


64 


P4, A-1 


K+ 


64.54 


2.569 


4282.4 


2.569 




bs 


P3, A-i 


bi 


Pi, A-i 


F- 


-11.45 


0.4304 


63.95 


1.042 



where 



*d(di,d2) : 



"^12 



di ■ d2 



3(di ■ri2)(d2 ■ r2i) 



r2 



and [1 — K{roi/a)] is the screening function [211231 . 

The repulsion between the electron shells of the ion and oxygen 
of the water molecule is modeled by the exponential functions (2] 






no 



J = l,2 



nj 



(D4) 



where 63 stands for the amplitude of the "ion-oxygen" repulsion 
energy, p3 is the inverse radius of action of the repulsion forces 
between the shells of the ion and oxygen, 64 is the amplitude of 
the "hydrogen-ion" repulsion energy, and p4 is the inverse radius of 
action of the repulsion forces between hydrogens and ion electron 
shells. Their values are gathered in Table 9. 

At distances significantly exceeding the sizes of a water molecule 
and the ion, the interaction potential takes the asymptote 



* -5> <I>d(dw,di), 



(D5) 



where dw = djj -I- do is the dipole moment of the water molecule 
and djj is the contribution made to this moment by hydrogens. 

1. F.H. Stillinger and C.W. David, J. Chem. Phys. 69, 1473 
(1978). 

2. I.V. Zhyganiuk, Ukr. Fiz. Zh. 56, 225 (2011). 



900 



ISSN 2071-0191 Ukr. J. Phys. 2011. Vol. 56, No. 9 



SPECIFIC FEATURES OF MOTION OF CATIONS AND ANIONS 



3. J.M.G. Barthel, H. Krienke, and W. Kunz, Physical 
Chemistry of Electrolyte Solutions (SteinkopfF, Darm- 
stadt, 1998). 

4. L.A. Bulavin, T.V. Karmazina, V.V. Klepko, and 
V.I. Slisenko, Neutron Spectroscopy of Condensed Media 
(Akademperiodyka, Kyiv, 2005) (in Ukrainian). 

5. D. McCall and D. Douglass, J. Phys. Chem. 69, 2001 
(1965). 

6. G.I. Safford, P.S. Leung, A.W. Naumann, and P.O. Chaf- 
fer, J. Chem. Phys. 50, 4444 (1969). 

7. P.S. Leung and G.I. Safford, J. Phys. Chem. 74, 3696 
(1970). 

8. L.A. Bulavin, A. A. Vasylkevych, and A.K. Dorosh, Ukr. 
Fiz. Zh. 31, 1703 (1986). 

9. L.A. Bulavin, G.M. Verbins'ka, L.O. Komarova, and 
V.V. Krotenko, Ukr. Fiz. Zh. 50, 939 (2005). 

10. L.A. Bulavin, M.P. Malomuzh, and K.M. Pankratov, 
Dopov. Nats. Akad. Nauk Ukr. No. 8, 70 (2009). 

11. A. Lyubartsev and A. Laaksonen, J. Phys. Chem. 100, 
16410 (1996). 

12. A. Lyubartsev, K. Laasonen, and A. Laaksonen, J. Chem. 
Phys. 114, 3120 (2001). 

13. A.V. Egorov, A.V. Komolkin, V.I. Chizhik, P.V. Yush- 
manov, A. P. Lyubartsev, and A. Laaksonen, J. Phys. 
Chem. B 107, 3234 (2003). 

14. O.Ya. Samoilov, Structure of Water Electrolyte Solutions 
and Hydration of Ions (AN SSSR, Moscow, 1957) (in 
Russian) . 

15. I.R. Yukhnovskii and LI. Kurylyak, Electrolytes 
(Naukova Dumka, Kyiv, 1988) (in Russian). 

16. M.F. Golovko, in Binary Distribution Functions for Sys- 
tems of Interacting Charged Particles (Evrosvit, Lviv, 
2010) (in Ukrainian). 

17. P. Blanckenhagen, Ber. Bunsenges. Phys. Chem. 76, 891 
(1972). 

18. T.V. Lokotosh, N.P. Malomuzh, and K.N. Pankratov, 
J. Chem. Eng. Data 55, 2021 (2010). 

19. Y. Kataoka, H. Hamada, S. Nose, and T.J. Yamamoto, 
Chem. Phys. 77, 5699 (1982). 

20. T.V. Lokotosh, N.P. Malomuzh, and V.L. Zakharchenko, 
Zh. Strukt. Khim. 44, 1104 (2003). 

21. L.A. Bulavin, A.I. Fisenko, and N.P. Malomuzh, Chem. 
Phys. Lett. 453, 183 (2008). 

22. L.A. Bulavin, T.V. Lokotosh, and N.P. Malomuzh, 
J. Mol. Liq. 137, 1 (2008). 

23. I.V. Zhyganiuk, Dopov. Nats. Akad. Nauk Ukr. No. 8, 77 
(2009). 

24. E.R. Nightingale, J. Phys. Chem. 63, 1381 (1959). 

25. J.C. Koneshan, J.C. Rasaiah, R.M. Lynden-Bell, and 
S.H. Lee, J. Phys. Chem. B 102, 4193 (1998). 



26. J.E. House, Inorganic Chemistry (Academic Press, San 
Diego, 2008). 

27. S.H. Lee and J.C. Rasaiah, J. Chem. Phys. 82, 1239 
(1985). 

28. I.D. Zaitsev and G.G. Zaitsev, Physical- Chemical Proper- 
ties of Binary and Multicomponent Inorganic Substances 
(Khimiya, Moscow, 1988) (in Russian). 

29. V.Ya. Antonchenko, A.S. Davydov, and V.V. Il'in, Fun- 
damentals of the Physics of Water (Naukova Dumka, 
Kiev, 1991) (in Russian). 

30. D. Eisenberg and W. Kauzmann, The Structure and 
Properties of Water (Oxford University Press, New York, 
1969). 

31. M.N. Buslaeva and O.Y. Samoylov, in The Chemical 
Physics of Solvation. Part A, ed. by R. Dogonadze et al. 
(Elsevier, Amsterdam, 1985). 

32. I.V. Zhyganiuk, Dopov. Nats. Akad. Nauk Ukr. No. 11, 
72 (2009). 

33. Yu.V. Ergin, O.Ya. Koop, and A.M. Khrapko, Zh. Fiz. 
Khim. 54, 392 (1980). 

34. E. Clementi and H. Popkie, J. Chem. Phys. 57, 1077 
(1972). 

35. H. Kistenmacher, H. Popkie, and E. Clementi, J. Chem. 
Phys. 61, 799 (1974). 

36. J. Frenkel, Kinetic Theory of Liquids (Dover, New York, 
1955). 

37. L.A. Bulavin, Doct. Thesis (Kiev, 1991) (in Russian). 

Received 01.07.11. 
Translated from Ukrainian by H.G. Kalyuzhna 



OCOBJIHBOCTI Pyxy KATIOHIB I AHIOHIB 
B P03HHHAX EJIEKTPOJIITIB 

JI.A. ByAaeiH, LB. XmauioK, M.Lt. MaAOMyatc, 
K.M. naHKpamoe 

P e 3 lo M e 

^ocjii.z],>KeHO cJ^isHHHy npHpo^y pyxjiHBOCTi ioniB i MOjieKyji bo^h 
y po36aBjieHHx bo^hhx posMnnax ejieKTpojiiTiB, kojih na o^hh Ioh 
npnna^ae He Sijitme n'5iTHaflii,5iTH MOjieKyji BOflH. lloKasaHO, iii,o 
noBefliHKa Koe(J)iij,ieHTiB pyxjiHBOCTi MOJieKyji BOflH i ioHlB, a TaKO>K 
Koe4)iij,ieHTiB caMOflH4>y3ii MOJieKyji BOflH BHpimanbHiiM hhhom bh- 
3HaMaeTbC5i pafliycaMH ix TBepflHx o6ojiohok, a ne bojiubom cItkh 
BOflHeBHx 3B'5i3KiB y CHCTeMi. BcTaHOBjieHO, iii,o BnjiHB ri^paTai^iii- 
HHx ecJjeKTiB na 3HaMeHH5i rycTHHn CHCTeMH i Koe4)iii;ieHTiB caMO- 
^H4)y3ii MOJieKyji bo;i,h He nepeBHmye kIjibkox BiflCOTKiB. Ha ocho- 
Bi MiKpocKoniHHiix yHBjieHB noKa3aHO, iii,o Bi;i,MiHHa noBe^inKa Ka- 
Tiona K^ Ta anioHa F~, mo MaiOTB o;i,HaKOBi jKopcTKi pa^iycii, .zi,o- 
6pe y3rofl>KyeTBCH 3 oco6jihboct5imh Mi>KMOjieKyji5ipHOi B3aeMO^ii, 
5iKa onHcyeTBCH ysarajiBHeHHM noTeHLi,iajiOM CTijiiH;i,>Kepa— /l,eBifla 

mi 
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